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V 12 - 50^ 



Ordinary B*geel Functions (SADOI Only) 

Closed 

2k2 

0-10 inclusive 

At least 7 decimal places 

Approximately 170 in. sec. 

3Ms routine calculates the Ordinary Bess el Functions which 

are defined as follows? 



m=0 ^n+2m I. , , v I 
2 m: (m+n) * 



2 (z) cos n* - J (z) 
N (z) * -£ 2 



n 



sin nit 



wn|ere n is an integer and is confined to M = 0, 1 in this 
routine, Th*se functions Sfre sdlutions to Bessel's equation: 



z* 



d w 

7? 



+ Z 



%- + (2* - n 2 ) w - 



dz 



The accuracy of the routine depends upon the accuracy of 
the argument. For this reason, a double precision argument 
may be used. In cases where the argument, scaled by a 
suitaola scaling factor, has sufficient accuracy in single 
precision, a single precision argument may be used'../ 
The true value of the argument must be greater than .000005. 
Values smaller than this will give incorrect results. 

To use a double precision argument, where the argument is 

-M 

scaled by 2 , store the most significant part in location 

OF and the least significant part in the accumulator and 
enter with 
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pfl 



50 MF 
50 pF 

26 () 



where M > but M < 39. 

If a single precision argument is to be used, the 
accumulator must be clear upon entry. One form for 
entry may be 



NOTE: 



p+1 



50 MF 

51 pF 



26 () 



Again the M3P must be in location OF. 

111© routine will leave the following? information in th& 

locations shown: 

J (x) 



in location 7F 

8F 

•" " 8F 

MM 



9F 
10F 
10F 



l/l6N (x) 

xs ± (x) 

N x (x) 



if x < 3 
if x > 3 

if x < 3 
if x > 3 



Notice the two forms for N Q (x) and N (x) depending on the 
size of the argument. 

The aim of the double precision argument is to allow the 
use of a wide range of arguments, all of which may be 
scaled by the same scaling factor. This means, for example, 
that if a range of arguments ia from .00005 to 2 12 , they 
may all be scaled by 2 J . The routine rescales each value 
down to the minimum value needed for each. For small 
arguments scaled by large factors, a double precision 
argument is needed to insure the accuracy of the routine. 
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This routine contains the subroutines; (S3) Natural 
Logarithm, (Rl) Square Root, (T5) Sine-Cosine. These may 
be used by other parts of the program and may be entered 
in the usual manner since they are referred to in this 
routine as symbolic addresses . 

The methods used to obtain these functions may be found 
in "Mathematical Tables and Aids to Computation" 195^ * 
Vol. 8, pp. 240-2^1. 
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i 


00K(V12) 




o 


40 IF 
ifl 132L 




1 


K5 F 






42 121L 


Plant link 


2 


10 20F 


") 




42 132L 


3 


41 133L 
126 197L 




4 


L5 F 






00 IF 


Rescale 


5 


. 36 6l 


v 




26 9L 


1 Argument 


6 


40 F 
F5 133L 




7 


40 133L 
LO 132L 


1 1 


8 


36 9L 
26 4L 




9 


L5 F 
4o 5F 


*S Arg. x 2" M 


10 


L5 132L 
LO 133L 




11 


40 6f 






L3 6F 


r> 


12 


32 24L 
40 F 


< Arg < 1 


13 


F5 F 






36 26L 


1 < Arg < 2 


14 


40 F 
F5 F 




15 


36 16L 


2 < Arg < 4 




26 74L 


Arg > 4 

' — 



LOCATION 
1 


ORDER 


\ l6 


L5 5F 


■ 


LO 122L 


17 


36 I8L 




26 28L 


18 


1*0 F 




L3 F 


19 


36 20L 




26 7^L 


20 


L5 123L 




hO 7F 


21 


L5 12^L 




40 8F 


22 


L5 125L 




40 9F 


23 


L5 126L 




40 10F 


2k 


22 121L 




L5 5F 


25 


10 2F 




26 27L 


26 


L5 5F 




10 IF 


27 


kO 5F 




50 iklL 


28 


L5 5F 




50 28L 


29 


26 (S3) 




10 6F 


30 


75 130L 




00 2F 


31 


LO 127L 




kO 131L 


32 


50 128L 




75 5F 
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Test for size 
of Argument 

2 < Arg < 3 



Arg = 3 
Arg > 3 



Arg = 3 

Place value of 

functions in 

locations 



Prepare for 
small arg 



l/l6jU> x/2 



LOCATION 



ORDER 



33 
34 
35 
36 
37 
38 

39 
40 
4l 
42 

43 
44 

45 
46 

47 
48 
49 



00 2F 
40 F 
L5 27L 
42 6JL 
L5 98L 
46 70L 
26 63L 
L5 196L 
L4 140L 
00 2F 
40 7F 
L5 73L 
42 67L 
L5 101L 
46 70L 
26 63L 
50 49L 
L5 196L 
L4 147L 
50 190L 
40 8F 
50 131L 
75 129L 
40 3F 
50 3F 

75 7F 
L4 8f 
40 8F 
L5 66L 
42 67L 
L5 4lL 
46 70L 
26 63L 
L5 196L 
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Arg/3 



Preset for 
J Q (small) 



J Q (small) 



Preset for 
N Q (small) 



l/l6 N Q (small) 

Preset for 
J, (small) 



LOCATION 



ORDER 



50 
51 
52 
53 
5* 
55 
56 
57 
58 

59 
60 
61 
62 
63 

a 

65 
66 

67 



Ik 15^L 
40 2F 
50 2F 
75 5F 
00 2F 
40 9F 
L5 62L 
1*2 67L 
L5 107L 
46 70L 
26 63L 
L5 196L 
Lfc 161L 
kO 10F 
50 88L 
50 3F 
75 9F 
ho 3F 
50 3F 
75 5F 
00 2F 
IA 10F 
00 k-F 
kO 10F 
22 121L 
50 162L 
50 F 
7J F 
kO IF 
kO 2F 
4l ;96L 
L5 135L 

40 136L 

50 155L 
50 2F 
75 F 
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J (small) 

Preset for 
^ (small) 



XN x (x) small 
Return to master 



Form sum 



LOCATION 



ORDER 



68 

69 
70 
71 
72 
75 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 



L4 196L 
40 196L 
F5 156L 
40 156L 
32 F 
50 IF 
75 2F 
40 2F 
F5 67L 
42 67L 
26 67L 
50 148L 
L5 6f 
42 79L 
42 95L 
42 99L 
L5 139L 
46 111L 
L5 138L 
42 90L 
50 134L 
L5 183L 
L5 122L 
10 F 
66 5F 
00 2F 
S5 F 
40 3F 
75 12 8L 
50 82L 
26 (Rl) 
40 4f 
L5 87L 
42 67L 



1 
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of terms 



Prepare for 
large Arg. 



i/YT 



LOCATION | 


ORDER 


NOTES PAGE 6 V 12 


85 


L5 57L 
46 70L 




86 


L5 3F 






40 IF 




8? 


40 2F 
50 l69L 




88 


26 65L 
50 176L 




89 


L5 196L 
L4 168L 




90 


40 10F 






L5 F 


Prepare for 


91 


42 67L 


... .■ ■-...,. ,.. .. : : . 




L5 117L 


large arguments 


92 


46 70L 






L5 3F 


, 


93 


40 IF 
40 2F 




9* 


26 65L 
L5 196L 




95 


L4 175L 
10 F 




96 


40 F 
L5 5F 




97 


LO F 
40 F 




98 


50 36L 
50 F 




99 


75 137L 
00 F 


■ 


100 


4o 9F 
iA 154L 




101 


50 4lL 






50 101L 


— : 



LOCATION 



ORDER 



102 

103 

104 

105 

106 

107 

108 

109 

110 

111 

112 

113 
114 

115 
116 

117 
118 



26 (T5) 

4o F 

50 10F 

75 4f 

40 6F 

50 6P 

75 F 

00 IF 

40 7F 

L5 9F 

50 55L 

50 H07L 

26 (T5) 

40 F 

50 F 

75 6f 

00 IF 

4o 8f 
22 F 

L5 138L 
46 111L 
L5 139L 
42 90L 
L5 78L 
42 67L 
L5 7F 
40 132L 
L5 8F 
40 133L 
26 85L 
50 94L 
LI 7F 
40 10F 
L5 8F 
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cos (Arg - 0) 



J Q (large) 
or J. (large) 



(T5) 

sin (Arg - 0) 



N Q (large) 
or N, (large) 



LOCATION 



ORDER 



119 


4o 9F 




L5 132L 


120 


ko 7F 




L5 133L 


121 


ko 8F 




22 F 


122 


00 F 




00 7500 0000 


123 


80 F 




00 7399 kdOk 


12^ 


00 F 




oo 3768 5000 


125 


00 F 




00 3390 5895 


126 


00 F 




00 32^6 7^0 


127 


80 F 




00 9566 7830 


128 


00 F 




00 3333 3333 


129 


00 F 




00 6366 1977 


150 


00 F 




00 6931 4718 


131 


00 F 




00 F 


132 


00 F 




00 F 


133 


00 F 




00 F 


13^ 


00 F 




00 F 


135 


LL 1+095F 




LL 4090F 



0000 J 

5100 J 

0000 J 

8530. J 

0000 J 
1220 J 

3333 J . 
2368 J 
O56O J 
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Return to master 



LOCATION ORDER 

■■•MWMMMlHi ' ; |- III' .. 

156 00 F 
00 F 

137 00 F 

00 3183 0988 6l84 J 

138 00 117L 
00 88L 

139 00 111L 
00 42L 

140 00 F 

00 2500 0000 0000 J 

141 80 F 

00 4375 0007 5000 J 

142 00 F 

00 3164 0520 0000 J 

143 80 F 

00 9209 0335 0000 J 

144 00 F 

00 111 1975 0000 J 

145 80 F 

00 9990 1390 0000 J 

146 OOF 

00 5250 0000 J 

147 00 F 

00 2296 6681 870 J 

148 00 F 

00 373 4960 3750 j 

. 149 80 F 

00 9535 3101 0000 J 

150 00 F 

00 158 1257 3120 J 

151 80 F 

00 9973 3675 1000 J 

152 00 F 

00 2 6744 75OO J 



NOTES 



PAGE 9 



V 12 



Coefficients 



for 



J (small) 



Scaled by 



,-2 



Coefficients 



for 



N Q (small) 



[scaled by 2 



-4 



LOCATION 


ORDER 
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153 


80 F 

Oo 9998 kkji 2500 J ^ 


J 


15^ 


00 F > 
00 5000 0000 0000 J 




155 


80 F 






00 4375 0015 0000 J 


Coefficients 


156 


00 F 






00 2109 3573 0000 J 


for 
> J (small) 


157 


80 F 




00 9604 5711 0000 J 




158 


00 F 






00 kk 3319 0000 J 


unsealed 


159 


80 F 






00 9996 8239 0000 J 


! 
i 


160 


00 F 

00 1109 0000 J i 


1 

i 




f 


161 


80 F 

00 9602 1126 4500 J 




162 


00 F 






00 138 2556 8750 J 


Coefficients 


163 


00 F 


for 




00 1355 I693 1200 J 


164 


80 F 






00 9177 1983 1300 J 


1 TH (small) 


165 


00 F 


1 scaled by 




00 195 2^69 3750 J 


2 A 


166 


80 F 

00 997J4- 9390 0000 J 




167 


00 F 






00 1 7^20 6250 J 


) 


168 


00 F 

' 00 7978 8^56 0000 J 


~\ 


169 


80 F 




i 


00 9999 9923 0000 J 


.Coefficients 



LOCATION 


ORDER 
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170 


80 F 


——————— — ' — _ — 1 




00 9944 7260 0000 J 


for 


171 


80 F 

00 9999 0488 0000 J 


- f o 


172 


00. F 






00 13 7237 0000 J 


unsealed 


173 


80 F 






00 9992 7195 0000 J 




17^ 


00 F 






00 14476 0000 J J 


175 


00 F 

00 7853 9816 0000 J 


1 


176 


00 F 






00 4i6 6397 0000 J 


Coefficients 


177 


00 F 






00 3954 0000 J 


for 


178 


80 F 


-*o 




00 9973 7^27 0000 J 


unsealed 


179 


00 F 

00 5 4l25 0000 J 




180 


00 F 

00 2 9333 0000 J 




181 


80 F 


: " 




00 9998 6442 0000 J 


J 


182 


00 F 

00 7978 8^56 0000 J 


^ 


183 


00 F 






00 156 0000 J 


Coefficients 


184 


00 F 

00 165 9667 0000 J 


I for 


185 


00 F 


f i 




00 1 7105 0000 J 


1 unsealed 


186 


80 F 

00 9975 0489 0000 J 





LOCATION 


ORDER 
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187 


00 F 

00 11 3653 0000 J 


J 






188 


80 F 

00 9997 9670 0000 J 


^ 






189 


00 F 

00 7853 9816 0000 J 








190 


80 F 

00 8750 0388 0000 J 


/ 






191 


80 F 

00 9999 ^+350 0000 J 


J Coefficients 






192 


00 F 

00 63 7879 0000 J 


1 for 






193 


80 F 

00 9992 5652 0000 J 




unsealed 






19^ 


80 F 

00 9992 0176 0000 J 










195 


00 F 












00 2 9166 0000 J 


J 






196 


00 F 
00 F 








197 


L3 132L 
36 9L 








198 


50 IF 
26 kL 

00K(S3) 


199-212 "~"\ 


Add these tapes 


in 




OOK(Rl) 


213-221 r 


this order. 






00K(T5) 


222-2J+2 / 








■■ ___j 


— . 







